Sign-time distribution for a random walker with a drifting boundary 
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We present a derivation of the exact sign-time distribution for a random walker in the presence 
of a boundary moving with constant velocity. 

There has been tremendous interest within the statistical physics community over the past five years in the per- 
sistence properties of various simple systems Q. Aside from the persistence probability itself, another quantity of 
great interest is the "sign-time" distribution (STD) [|||] (also known as the "occupation-time" or "residence-time" 
distribution). It seems that STD's are very difficult to calculate, even for the simplest systems. Famous examples for 
which the STD is known exactly are the pure random walk, yielding the arcsine law of Levy |^J5|, and the renewal 
process, originally studied by Lamperti Q, and recently revisited There have been concerted efforts recently 

to find the STD for a generalized random walk model ||[l0]], which have met with some success, although an exact 
solution is still not known. Here we consider the STD for a random walker in the presence of a boundary moving 
with constant velocity (which is equivalent to a random walker with a constant drift velocity in the presence of a 
stationary boundary). The first passage probability jflj and the persistence probability []12j for this problem are both 
known. However, to the author's knowledge, the STD has not been previously calculated. The purpose of this paper 
is to present a relatively straightforward derivation and analysis of this latter quantity. 

We consider a random walker whose position is denoted by x(t), which satisfies the equation of motion 

i = £(*), (1) 

with x(0) = 0, and (£(<)) = and (£(£)£(£')) = DS(t — t'). In addition we consider a moving boundary whose position 
is denoted by b(t). In the main part of this paper we shall consider linear motion of the boundary, i.e. b(t) = vt. 

The sign-time r^(t) is defined as the proportion of elapsed time t the walker has spent on the positive side of the 
boundary: 

t 

r i (t) = Jdt'e(x(t')-b(t')), (2) 



where Q(x) is the step function. We are interested in calculating the STD P(r,t) defined by 

P(r,t) = {S(r-r i (t))) . (3) 

By making an integral representation of P(r,t), we can reduce the problem to a calculation of the moments of the 
STD: 

OO 

— OO 
OO 

n=0 



— OO 



To this end we introduce a convenient double-integral representation of the step function 



OO OO 

dk 



Q(x) = I du5(x-u) = I du / — e lk{x - u) . (6) 



o o 



Utilizing this in the definition of the sign-time, we can write the n moment as the following time-ordered expression: 
<T 5 (tn =n\ I ■Dt(n)(Q{x(t 1 )-b(t 1 ))---Q(x(t n )-b(t n )) ) , (7) 
= n\ [ Vt{n) [ Vu(n) f Vk{n) e -<M«i+»(*0) . e -*M"»+*(*»)) / e ^*(*i) . . . e «*»*(*»)^ , (g) 
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where the multiple integral symbols are defined via 

t ti t„ 

Vt{n) = / dti ( di 2 ■ ■ ■ / dt n , (9) 







T>u(n) EE / dtti " ' / ^"n i (10) 


oo 

dki f dk n 



J Vk{n) 



2. / 2. • ^ 



Writing the solution of Eq.(|]) as x(t) = J dt'£(t'), the average on the right-hand-side of Eq.(|8|) is easily performed to 

o 

give 



^ e ifcx»(tO . . . e ik nX (t n )^ = exp |_| [ fc 2 (ti _ ta) + (fci + fe) 2 (<2 _ fa) + . . . + (fci + . . . + kn) 2 tn } | 



(12) 



Inserting this expression into Eq.(||) and changing variables via (pi — ki, . . . ,p n — ki + ■ ■ ■ + k n ), we find the integrals 
over {pi} decompose and may be performed, giving us 



(r c (t) n )=n! J Vt(n) J Vu{n) G{m + 6(t x ) - u 2 - b(t 2 ),h - t 2 ) ■ ■ ■ G(u n + b{t n ),t n ) , (13) 

where G(u, t) is the random walk propagator 

G(u,t) = (2nDt)~ 1/2 exp[-u 2 /2Dt] . (14) 

We shall now focus on the case of a boundary drifting with constant velocity v, i.e. we write b(t) — vt. We note the 
appearance of a fundamental rate constant 7 = v 2 /2D. It is convenient to introduce the function H (u, t) = G(u+vt, t). 
By considering the Laplace transform of the n th moment, and invoking the Laplace transform convolution theorem 
iteratively, we find 



Ti 

A>|t [(T^(t) n )} = — I T>u(n) H(ux - U2, s) ■ ■ ■ H(un-i - u n , s)H(u n , s) , (15) 



where the Laplace transform operator is denoted by 

f( S )=£ s M(t)}= I dte- st f{t) . (16) 



Now, H(u, t) = cxp(—uv/D — jt)G(u, t) and consequently H(u, s) = e wv / D G{u, s + 7). The Laplace transform of 
G(u,t) is given by 

G(u,s) = (2sD)- 1 ^cxp[-(2s/D) 1 / 2 \u\} . (17) 
Substituting the explicit form for H (u, s) into Eq.(|T5|) and rescaling the integration variables we find 

£e\t [fo(*) w )] = s r 2(s t 7 )l" J V < n ) e ' a{s)Ul exp(-K- U2 | K_x - u n \ - \u n \) , (18) 

where we have defined a(s) = ( 7 /(s + 7 )) 1 ^ 2 - Eq. (|l8l) can be conveniently rewritten, for n > 1, as 



£*\t [fo(*) n >] = a[2 ( s + 7)]" ' dui e ' a{s)ulJ ^M > ( 19 ) 



2 



where the function J n (u) satisfies the integro- difference equation 

oo 

J n (u) = fdu' e-\ u - u '\j n ^(u') 



(20) 



with Jo = e u . This equation for J n (u) can be solved using generating function methods, and the details are relegated 
to an Appendix. The end result is 



f dz (\- (1- 2z) 1 / 2 \ r , N1 n 
Jn{u) = f^n { z"+2 ^ j exp[-(l - 2*)^] , 
c 



(21) 



where the contour C encircles the origin. 

We now substitute this solution for J n (u) into Eq.([l9]) and perform the integral over u\. Using the fact that 
§ c dz/z n+1 = for n > 1 allows us to rewrite the integrand so as to obtain 



£s\t[(n(t) n )} = 



n!(l + a(s)) f dz 



[2(s + 7)] rl J 2m z n+1 [a(s) + (1 - 2z) 1 /2] 



(22) 



The integral around the contour C may be re-expressed as an integral across the cut emanating from the branch point 
at z = 1/2 with the result 



C\t [(rdt) n )} = 



n!(l + a(s)) 



dq 



.1/2 



7ts(s + 7)" 7 (1 + <7)' i+1 (a(s) 2 + g) 



(23) 



It can be explicitly checked that this expression is also correct for n = 0, namely, that it gives the result 1/s. 

It is possible to invert the Laplace transform at this point, but it is more convenient to postpone this operation 
and instead introduce the integral representation for n! via 



nl = dy y n e~ v 



(24) 



Then Eq.(E3[) takes the form 



Referring to Eq.Q we have 



(! + «(*)) 



ITS 



dq 



a/2 



(l + q)(a{s) 2 + q) 



dy e v 



(a + 7)(l + <7) 



2tt 



IU) ) 



(25) 



(26) 



Substituting Eq.(p5|) into Eq.([26|), we can explicitly perform the sum over n which allows us to complete the integral 
over uj yielding a Dirac <5-function with argument (r — y/(s + "/)(! + q)), which in turn allows us to trivially perform 
the integral over y. This leaves us with (after rescaling the q integration variable), 



C a \ t [P(j,t)] 



OO 

= a( a )(l + q( a ))( 8 + 7) r _ f ^ T g^e-™ 



+ 



(27) 



It is worth noting the clean disappearance of the Laplace transform variable s from inside the g-integral, which allows 
us to complete the calculation easily from this point. Formally inverting the Laplace transform gives us 



a(a)(l + a(fl))(fl + 7) c - aT 



(28) 
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where we have defined 



N(0) EE dq 

TT 



qV 2 e-Pi 



,-0 



erfc(/3 1/2 ) , 



(29) 



(1 + q) (7T/3)V2 
o 

with erfc(z) the complementary error function fll3| ]. On performing the inverse Laplace transform in Eq.(p8|) we obtain 

P(T,t)= 1 N( 1 T)[2 + N( 1 (t-r))] . (30) 

This may be cast into a more symmetrical form by defining the function F + and its companion F~ via 

F ± (P) = T (7r/3) 1 / 2 erfc(±/3 1/2 ) . (31) 

We then have our final result in the form 

F+( 1T )F-(j(t-T)) 



P(r,t) = 



n(T(t - t))V2 



(32) 



It is interesting to note that the distribution is a product of a function of r and a function of (t—r). It is often convenient 
to express the STD in dimcnsionless variables. Defining </> = r/t G (0,1) and 77 = we have P(t, t)dr — ^(<p,rj)d(j) 
and consequently 



ri) 



F+(^)F-(»;(1-^)) 



7r(^(l-^))Va ^ 

We briefly examine some limits of the STD. For small velocity, or small times, we have 77 <C 1, and thus 



7^(1 - 



1/2 



1 



1 



(1 - 0)Va ^i/s 



1 



r(0(l - £))V2 



- 1 



0(?? 3/2 ) 



(33) 



(34) 



where the leading term is the arcsine law of Levy, as required. 

The limiting forms of W for 0-^0 and (f> — > 1 are given below, along with the various forms of large-77 behaviour: 



ttS 1 / 2 ' 



F+(ri) 



(A) 



1/2 



2tt7 ) (1-0) 1 /2 



0^0 , 

0-1, 

77 — > 00, rj<j) — » cx), 77(1 — 0) — > 

77 — > 00, 0^0, 77^ — > , 

rj — * 00, — ► 1, 77(1 — </>) — > 



(35) 



The distribution has a single minimum whose position </> m i n (?7) approaches unity as </> m in(?7) ~ 1 — r/77 for ij > 1, 
where the universal number r satisfies the transcendental equation 



( 7 rr) 1 / 2 e r erfc(-r 1 / 2 ) = (1 - 2r)/2r , 



(36) 



and has the value r = 0.2040539 . . .. 

It is hoped that the method of derivation presented here may be of use in calculating STD's for other simple 
stochastic processes, and that the explicit solution given in Eq.(|2]) will find useful application in a range of problems. 
It is currently being employed in an ecological context, to calculate survival probabilities for organisms in the presence 
of a step-like drifting environmental boundary JlJ]. 

The author would like to thank Satya Majumdar for an enjoyable correspondence, and for generously sharing his 
expert knowledge of the literature. 
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APPENDIX A: 



In this appendix we calculate the function J n (u), which satisfies the integro-difference equation 

oo 

J n {u) = fdu' e- |W| J„-i( M ') , (Al) 



with Jq(u) — e u . This is accomplished by means of the generating function 

oo 

J(u,z) = J2 znj n(u) • 



(A2) 



n=0 



Summing over Eq.(Al) we find 



J(u, z) = e~ u + z du' e -l W l J(u', z) . 



Differentiating twice with respect to u yields the differential equation 

dlJ(u,z) = (I - 2z)J(u,z) 
Demanding that J{u, z) is finite as u — > oo gives 

•f-(l-2z) 1 /2 



J{u, z) 



^exp[-(I-2z) 1 / 2 i 



(A3) 



(A4) 



(A5) 



where the z-dependent prefactor is found from substitution into Eq.(A3). The function J n (u) may be recovered from 
the generating function by means of the contour integral 



J n {u) 



dz J{u,z) 
2tt7 z n+1 ' 



(A6) 



where the circular contour C runs counter-clockwise about the origin and excludes all singularities other than the 
pole at the origin. 
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